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ABSTRACT 

We present an analysis of the star formation history (SFH) of the transition-type 
(dIrr/dSph) Local Group galaxy LGS-3 (Pisces) based on deep photometry obtained 
with the Advanced Camera for Surveys onboard the Hubble Space Telescope. Our obser- 
vations reach the oldest main sequence turn-offs at high S/N, allowing a time resolution 
at the oldest ages of a ~ 1.1 Gyr. Our analysis, based on three different SFH codes, 
shows that the SFH of LGS-3 is dominated by a main episode ~ 11.7 Gyr ago with a 
duration of ~ 1.4 Gyr. Subsequently, LGS-3 continued forming stars until the present, 
although at a much lower rate. Roughly 90% of the stars in LGS-3 were formed in 
the initial episode of star formation. Extensive tests of self-consistency, uniqueness. 
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and stability of the solution have been performed together with the lAC-star/IAC- 
pop/MinnlAC codes and these results are found to be independent of the photometric 
reduction package, the stellar evolution library, and the SFH recovery method. There is 
little evidence of chemical enrichment during the initial episode of star formation, after 
which the metallicity increased more steeply reaching a present day value of Z ~ 0.0025. 
This suggests a scenario in which LGS-3 first formed stars mainly from infalling fresh 
gas, and after about 9 Gyr ago, from a larger fraction of recycled gas. The lack of 
early chemical enrichment is in contrast to that observed in the isolated dSph galaxies 
of comparable luminosity, implying that the dSphs were more massive and subjected to 
more tidal stripping. 

We compare the SFH of LGS-3 with expectations from cosmological models. Most 
or all the star formation was produced in LGS-3 after the reionization epoch, assumed 
to be completed at z ~ 6 or ~ 12.7 Gyr ago. The total mass of the galaxy is estimated to 
be between 2 and 4 x 10® Mq, corresponding to circular velocities between 28 km s~^ to 
36 km s~^. These values are close to but somewhat above the limit of 30 km s~^ below 
which the UV background is expected to prevent any star formation after reionization. 
Feedback from SNe associated with the initial episode of star formation (mechanical 
luminosity from SNe = 5.3 x 10^* erg s^^) is probably inadequate to completely 
blow away the gas. However, the combined effects of SN feedback and UV background 
heating might be expected to completely halt star formation at the reionization epoch 
for the low mass of LGS-3; this suggests that self-shielding is important to the early 
evolution of galaxies in this mass range. 

Subject headings: galaxies: dwarf, galaxies revolution, galaxies: photometry, galaxies: stellar 
content, galaxies:structure, cosmology: early universe 



1. Introduction 

The present paper is part of the Local Cosmology from Isolated Dwarfs (Lcrojll project. 
The Hubble Space Telescope (HST) has been used to obtain deep photometry of six isolated dwarf 
galaxies in the Local Group: IC 1613, Leo A, Cetus, Tucana, LGS-3 and Phoenix. Color-magnitude 
diagrams (CM P) reaching the o ldest turn-off points have been derived. The Advanced Camera for 



Surveys (ACS. lFord et al.lll998l ) has been used to observe the first five, while Phoenix was observed 
with the Wide Field and Planetary Camera-2 (WFPC2). The main goal of the LCID project is 
to derive the star formation history (SFH), age- metallicity relation (AMR), and stellar population 
gradients of this sample of galaxies. Our objective is to study their evolution at early epochs and to 
probe effects of cosmological processes, such as the cosmic UV background subsequent to the onset 



^Local Cosmology from Isolated Dwarfs: [http: / /www.iac.^ 
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of star formation in the universe or physical processes such as the gas removal by supernovae (SNe) 
feedback. Our sample consists of field dwarfs which were chosen in an effort to study systems as 
free as possible from environmental effects due to strong interactions with a host, massive galaxy. 
The SFH is a powerful tool to derive fundamental properties of these systems and their evolution, 
but to study the earliest epochs of star formation, deep CMDs, reaching the oldest main sequence 
turn-offs, are required. 

Here we present our analysis of LGS-3. LGS-3 is a dIrr/dSph low-luminosity galaxy discovered 
bv iKarachentseval (119761). From its locati on, it appears to be a member of the Andromeda subgroup 
of the Lo cal Group (Ivan den Berghll200Cll). The dis t ance to the ga laxy has been estimated by several 



authors. iLed (|l995l ). lTikhonov fc Makaroval (|l996l ). lMould (|l997l ). and lAparicio. Gallart. fc Bertelh 



(jl997bl ) used the tip of the red giant branch (TRGB) to estimate distances ranging from 0.6 to 
0.96 Mpc. This me asurement is i nhere ntly challenging because the TRGB region of LGS-3 is very 
sparsely populated. iMiller et al.l (j200ll ) also used the TRGB but complemented with observations 
of the horizontal branch (HB), the red clump (RC), and the overall match of the CMD to model 
CMDs to estimate a d istance of 0.62 it 0.02 Mpc. Finally, from a study of the RR-Lyr stars in 
LGS-3, iBernardI bood ) derived a distance of 0.65 ± 0.05 Mpc ((m - M)o = 24.07 ± 0.15 mag). We 
will assume this value in this paper. 

Th e stel lar content o f LGS - 3 has b e en st udied by Lee ( 1995 ). Mouldl ( 1997 ). Aparicio et al 



(|l997bl ). and lMiller et al. 



(2001 



Mouldl (119971 ) discovered the first Cepheid variable candidates in 



this galaxy. These works show that LGS-3 contains stars of all ages; it has been forming star s since 
at least 12 Gyr ago and shows a low re cent star formatio n rate (SFR) (lAparicio et al.lll997bl ). The 



SFH shows a dependency with radius (jMiller et al 



2001 



)• 



No H II regions have been found in the galaxy ([Hodge &: Milleijll995l ). which is consistent with 
the low recent star formation activity. The m etalli c ity of the stars has been estimated from the color 
of the red giant bra nch (RGB). In t his wa y, iLed (jl995l ) obtained a mean metallicity of [Fe/H] = 
-2.10 ± 0.2 2 while l Aparicio et ahl (jl997bl ) derived a metallicity range -1.45 < [M/H] < -1.0. 

Miller et al.l koQll ) derived [Fe/H] = -1.5 ± 0.3 for the stars older than 8 Gyr, and [Fe/H] 1 

for the most recent generation of stars. 



The HI content of LGS-3 was first studi ed with single dish s pectra bv lThuan Sz Martini (jl979l ) 
and, more recently, with interferometry by lYoung fc Lol (| 19971). Using the result s by the latter 
authors and the adopted distance of 0.65 Mpc, the results from lYoung Sz Lol (|l997l ) correspond to 
a total HI mass of Mhi = 2.7 x 10^ Mq, which corresponds to a total atomic gas mass of 3.8 x 10^ 
M.Q when corrected for helium. They also measured a gas velocity dispersion of 8 km s^^ and 
negligible rotation at a galactocentric radius of 470 pc. The corresponding virial mass estimate is 

Myir = 2.1 X 10"^ Mq. 

In this paper, we present the SFH of LGS-3 obtained from observations with the ACS on the 
HST. The photometry reaches the oldest main sequence turn-offs of the galaxy, allowing us to obtain 
an accurate SFH even for the oldest stellar populations. The structure of the paper is as follows: in 
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^the observations and data reduction are discussed, the CMD is presented, and the completeness 
tests are described. In ^|3]we describe many of the details of the lAC-pop/MinnlAC methodology 
and its application to the LCID program observations. The derived SFH of LGS-3 is presented 
in ^ and is compared with those of other LCID galaxies in ^ The consequences of SNe driven 
winds and cosmic reionization on the SFH of LGS-3 are discussed in f|6l The main conclusions of 
the work are summarized in f|3 Finally, in appendix El the uniqueness, self-consistency, stability, 
and statistical significance of the solution are discussed. 



Observations and Data Reduction 



The ACS observations of LGS-3 were obtained between September 12 and 17, 2005. The 
F475W and F814W bands were selected as the most efficient combination to trace age differences 
at old ages, since they provide the smallest relative error in age and metallicity in the main-sequence 
and sub-giant regions. Total integration times were 15072 s in F475W and 13824 s in F814W. The 
observations were organized in 6 visits of 2 orbits each, and each orbit was split into one F475W 
and one F814W exposure. In the first orbit of each visit, exposure times were 1250 and 1147 s 
respectively while in the second orbit they were 1262 and 1157 s respectively. Dithers of a few pixels 
between exposures were introduced to minimize the impact of pixel-to-pixel sensitivity variations 
("hot pixels") in the CCDs. The observed field o f LGS-3 is sh own in Figure [H The core radius of 
LGS-3 is 0;82 and its optical scale length is 0'.78 (jMateolllQOsI ). so the 3'A x 3;4 format of the ACS 
covers a substantial fraction of the stars in LGS-3. Table [J gives a summary of the observations. 

We analyzed the images taken directly from the STScI pipeline (bias, fiat-field, an d image 
disto rtion corrected). Two PSF-fittiii g photometry packages, DAOPHOT/ALLFRAME (Istetson 
19941 ) and DOLPHOT (iDolphinI l2000l ) . were used independently to obtain the photometry of the 
resolved stars. Non-stellar objects and stars with discrepant large uncertainties were rejected based 
on estimations of profile sharpness an d goodness of fit. The final photometric lists contain ~ 30, 000 

(j2010bl ) for more details about the photometry reduction procedures. 



stars each. See Monelli et al. 



Individual photometry catalogs were calibrated using the equations provided by lSirianni et al 

(|2005l ). The zero-point differences between the two sets of photometry are ~ 0.05 for F475W- 
F814W and ~ 0.03 for F814W. These small differences are ty pical for obtaining HST photometry 
with different methods (iHill et al.lll998l : iHoltzman et al.ll2006l ). In fj3]we discuss a technique which 



minimizes the already small impact of these differences in the obtained SFHs. 

The CMD of LGS-3 is shown in Figure [2] for the DOLPHOT photometry. Individual stars 
are plotted in the left panel and density levels are shown in the right panel. The left axis shows 
magnitudes in the ACS photometric system corrected for extinction. Absolute magnitudes are 
given on the right axis using the adopted values for the distance modu lus ((m — M)o = 24.07) and 
extinctions {ApAj^w = 0.156 and ^_f814w) = 0.079 mag) (lBernardll2009l ). In the left panel, the lines 
across the bottom show the 0.25, 0.50, 0.75, and 0.90 completeness levels derived from false star 
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tests (discussed later in this section). Finally, in order to highlight the main featu res of the CMD 



three isochrones from the BaSTI stellar evolution library (jPietrinferni et al.ll2004l ) are also shown 
in the right panel. 

The CMD of LGS-3 shown in Figure [2] reaches ~ 1.5 mag below the oldest main-sequence 
turn-off : these observation s are ~ 2 mag fainter than the deepest CMD previously obtained for this 
galaxy (|Miller et ahlboOlh . A comparison of the observations with the over-plotted isochrones and 



the presence of an extended HB indicate that a very old, very low metallicity stellar population is 
present in the galaxy. The gap produced in the HB by the RR-Lyrae variables is visible as well as 
a narrow RGB. The latter is populated by stars all older than 1-2 Gyr and its narrowness is an 
indication of a low spread in metallicity. The RC can be observed at the red end of the HB. A MS 
with stars younger than 200 Myr is clearly apparent. Some blue and red core helium-burning stars 
might be also present above the RC {msu < 22.3; 0.5 < (777,475 — ttisu) < 1.2). The RGB bump can 
be also observed at (777,475 ~ "7814) ~ 1.4 and 777814 ~ 23.2 mag. This evolutionary f eature may be 



used a s indicator of a dominant coeval stellar population homogeneous in metallicity (IMonelli et al 



2010a|). 



Signal-to-noise limitations, detector defects, and stellar crowding can all impact the quality of 
the photometry of resolved stars with the resulting loss of stars, changes in measured stellar colors 
and magnitudes, and systematic uncertainties. To characterize these observational effects, we have 
used the standard technique of injecting a list of false stars in the observed images and obtaining 
their photometry in an identical manner as for the real stars. These observational effects must be 
simulated in the stars of the global, synthetic CMD (sCMD) (hereafter synthetic stars) to be used 
for the analysis of the SFH (see below). 



Monelli et al.l (j2010bl ) provide a detailed description of the procedures we adopt for the char- 
acterization and simulation of these observational effects. For the present case of LGS-3, ~ 5 x 10^ 
false stars were used for the DOLPHOT photometry and ~ 6.5 x 10^ were used for the DAOPHOT 
photometry. The completeness factor, F, is defined as the rate of recovered to injected stars as 
a function of magnitude and color. Figure [21 right panel, shows F = 0.25,0.50,0.75 and 0.90 for 
the DOLPHOT photometry. Results for DAOPHOT are similar. At the magnitude level of the 
oldest MS turn-off stars, F ~ 0.90, allowing for very strong constraints on the oldest epochs of star 
formation. 

We simulate the observational effects in the stars of the sCMD with the routine obsersin. It uses 
the unrecovered false stars and the differences between the injected (777,) and measured magnitudes 



(777^) of the recovered stars for the simulation (jAparicio fc Gallartlll995l ). Since the density of stars 



(and thus the amount of crowding) normally varies substantially across the galaxy's field, spatial 
information has also been taken into account in the completeness test. To do so, first the spatial 
density distribution of observed stars in LGS-3 is measured. For that, a grid of 100 x 100 pixel size 
bins is defined, stars are counted in each bin and the result is smoothed and normalized. This is 
used as a density probability distribution to assign positions to each synthetic star. Then, for each 
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synthetic star with magnitude m<j, a hst of the false stars with |mj — ms\ < Cm', r < Cr is created for 
each filter, where r is the spatial distance between synthetic and false star and €m and free input 
searching intervals. If the number of false stars selected in this way, in common to both filters, is 
< 10, the parameters €m and are increased in turn until 10 or more false stars are found. Then, 
one of them is randomly selected from this list. If that star was unrecovered, the synthetic star 
is eliminated from the sCMD. If the selected star was recovered and m[ and m'j. are its injected 
and recovered magnitudes in a given filter, then = + — will be the magnitude of the 
synthetic star with observational effects simulated. The same is done for both filters. 



3. The lAC Method to Solve for the SFH 

3.1. Details of the lAC-star, lAC-pop, and MinnlAC codes 

Our analysis of the SFH of LGS-3 and the rest of LCID galaxies is done upon two fundamental 
requirements: i) the consistency of the method, in order to obtain SFHs that are readily comparable 
and ii) the stability and uniqueness of the solution. To accomplish i), we have adopted the same 
criterion on the selection of input parameters to derive the SFH for all galaxies of the LCID project. 
For ii), we have further developed our method and created new tests and control procedures. Details 
are given in the following and appendix 1X1 



The lAC method to solve the SFH is based in three main codes: (i) lAC-star (jAparicio &: Gallart 



2004 ) is used for the computation of synthetic stellar populations and CMDs; (ii) lAC-pop (jAparicio &: Hidalgo 



20091) is the core algorithm for the computation of SFH solutions, and (iii) MinnlAC (introduced in 



Aparicio Sz Hidalgdl2009l . but explained in more detail in this paper) is a suite of routines created to 
manage the process of sampling the parameter space, creating input data, and averaging solutions. 
As a summary. Figure [3] shows a data fiow diagram for the SFH computation process. 

Considering that tiin e and metallicity are the most important variables in the problem (see 



Aparicio Sz Hidalgo! l2009l ) . we define the SFH as an explicit function of both. We will use two 
definitions, one in terms of number of formed stars and another in terms of gas mass converted 
into stars. Formally, ipNit, Z)dtdZ is the number of stars formed within the time interval [t, t + dt] 
and within the metallicity interval \Z, Z + dZ], while V'(i) Z)dtdZ is the mass converted into stars 
within the same intervals. The relation between ipNit, Z) and ip{t, Z) is the average stellar mass at 
birth. Both can be identified with the usual definitions of the SFR, but as a function of time and 
metallicity. 

Several criteria can be used to grid the CMD. The simplest is to use a uniform grid and count 
stars within each box. However, not all the regions of the CMD and stellar evolutionary phases 
contain information of the same quality for the purpose of the SFH derivation. An alternative is to 
give higher weights to regions with more certain stellar evolutionary phases (like the main sequence 
or the sub-giant branch), regions which are more heavily populated, or regions in which the photom- 
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etry is more accurate. To do so, we have introduced what we call bundles (see lAparicio fc Hidalgo 



20091 . and Figure H]), that is, macro-regions on the CMD that can be subdivided into boxes using 
different appropriate samplings. The number and sizes of boxes in each bundle determines the 
weight of that region in the derived SFH. This is an efficient and flexible approach for two reasons. 
First, it allows a finer sampling in those regions of the CMD where the models are less affected by 
uncertainties in the input physics. Second, the box size can be increased, and the impact on the 
solution decreased, in those regions where the small number of observed stars could introduce noise 
due to small number statistics or where we are less confident in the stellar evolution predictions. 

lAC-pop provides a solution for the input set of data and parameters. It can also provide several 
solutions for the same set by changing the initial random number generator seed. Alternatively, 
it can deliver several solutions by randomly changing the input data (the number of stars in 
the original, oCMD, boxes) according to a Poisson statistics. However all of these solutions are 
eventually sensitive to several other input functions and parameters. To better analyze their effects, 
we can divide them into four groups. The first group consists of sampling parameters; the main 
ones are the number and boundaries of parent, pCMDs, and the criteria to define the bundles and 
grid arrays for sampling the CMDs. The second group consists of external parameters; these are 
the photometric zero points, distance modulus, and reddening. The third group consists of model 
functions, which are those functions defining the sCMD; the main ones are the initial mass function 
(IMF), (/>(m), and the function accouiiting f or the frequency and relative mass distribution of binary 



stars, f3{f,q) (see lAparicio &: Hidalgo! 12003 ). Finally, we have the external libraries, which consists 
of stellar evolution and bolometric correction libraries. Limiting and measuring the effects of all 
these parameters in the SFH solution requires several thousand runs of lAC-pop using different 
input sets. MinnlAC is used to test for the effects of the four groups are treated separately. In 
general, the solutions obtained with varying sampling parameters are averaged; their rms deviation 
is also a good estimate of the solution internal error. As for the external parameters and model 
functions, the set with the best Xu is usually adopted. Regarding the external libraries, the solution 
with the best xt may be adopted, but all the solutions are usually conserved; their differences can 
provide an indication of the external uncertainties associated to such libraries. 

MinnlAC (from Mmnesota and I ACT) is a suite of routines developed specifically to manage 
the process of selection of sampling parameters, creating input data sets, and averaging solutions. 
More in detail, regarding the parameter selection and data input, MinnlAC is used for two tasks. 
First, to divide the sCMD in the corresponding pCMD according to the selected age and metallicity 
bins. Second, to define the bundles and box grids in the CMDs and counting stars in them. With 
this information, MinnlAC creates the input files for each lAC-pop run. Several selections of both 
simple populations and bundle grids are used within each SFH solution process. MinnlAC does 
this task automatically, requiring only the input of bundles and corresponding grid sizes, age and 
metallicity bin sizes and the steps by which the starting points of the CMD grid as well as age 
and metallicity bins should be modified from one lAC-pop run to another. A second module of 
MinnlAC is used to average single solutions found by lAC-pop as required. 
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3.2. Solving procedure, results, and uncertainties for LGS-3 

In this section we will describe with some details the specific application of the lAC method 
to the case of LGS-3, detailing the different effects mentioned in ^ sampling parameters, external 
parameters, model functions and external libraries. Our choices are summarized in Table [2j For 
each item, Columns 2, 3, and 4 give minimum (or choice 1), maximum (or choice 2) and the step 
used (if relevant), respectively. Column 5 gives the final adopted value (if relevant). The SFH 
of LGS-3 has been derived independently for the DOLPHOT and DAOPHOT photometry sets 
following the same procedure as described below. 

Using lAC-star, several sCMDs were computed, one with each of the IMFs and binarity func- 
tions listed in Table [2j Each sCMD was built with 8 x 10^ stars with a constant star formation 
rate in the range of age 0-13.5 Gyr and a metallicity range 0.0001-0.005. The latter value has 
been assumed after a preliminary test with isochrones which showed that stars with metallicity 
larger than 0.005 do not exist in LGS-3. Using a large number of synthetic stars in the sCMD 
is important in order to minimize the stochastic effects of the Monte Carlo simulations used in 
lAC-star; to have a statistically significant number of stars in each pCMD (> 6000 stars), and to 
reduce the uncertainties in the stellar composition of simple populations in the pCMDs. 

The bundles and one of the grid sets used to bin the CMDs are shown in FigureHJ In total, 
72 combinations of IMF and binarity {model functions) have been tested. For each of them, 24 
solutions have been obtained by varying the CMD binning and the simple stellar population sam- 
pling (sampling parameters). The average -0 of these 24 solutions is adopted as the solution for 
the specific IMF and binarity. These averages have been obtained by using a boxcar of width 1 
Gyr and step 0.1 Gyr for t and width 0.0006 and step 0.0001 for Z. A Xu is also obtained as the 
average of the 24 single Xu values. The average solution with the smallest Xu is assumed as the best 
one. The IMF and the binary fraction corresponding to this solution are assumed here as the ones 
best representing our photometry data, although we do not intend to reach any strong conclusions 
concerning their true values. 

To minimize the effects of external parameters (distance, reddening and photometry zero- 
point), the SFH is derived for different offsets in color and magnitude applied to the oCMD. The 
offset giving the minimum Xu is assumed as the best. These are obtained when the oCMD are 
shifted by [A(F475W-F814W), A(F814W)] = [0.02,0.11] and [-0.03,0.11] for the DOLPHOT and 
DAOPHOT photometry sets, respectively. These tests also provide a way to estimate how distance, 
reddening and zero-point unc6rt8.inties translcite into the solution of the SFH. Defining X'^z^mm 
the xi value of the best solution, if a -2 is its corresponding rms deviation, all the solutions 

A, h',min 

with x^jy < X^umin + "^v'a should be considered valid. We can assume that the rms deviation 
of all the corresponding single solutions is a proxy of the error of the adopted solution, ap. 

The former error, ap, includes the effects of sampling and external parameters. However, 
there is one more source of uncertainty that must be taken into account, namely the effect of 
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statistical sampling in the oCMD. To estimate this, the solution is again obtained after varying the 
numbers of stars in the oCMD grid according to a Poisson statistic. This is done 20 times and the 
rms deviation of the solutions, cj/, is obtained. The final adopted error is the sum in quadrature 



Finally, solutio ns have been obtained using both the BaSTI (jPietrinferni et al.ll2004l ) and the 
Girardi et al.l (120001) stellar evolut ionary libraries {external libraries). In both cases the bolometric 
corrections by iBedin et al.l (j2005l ) were used. Both produced qualitatively similar results and a 
comparison between them together with the results obtained with other methods is presented in 
4.31 For detailed analysis presented in the following sections, for simplicity, we use only the BASTI 
stellar evolutionary libraries. 



4. The SFH of LGS-3 

4.1. Main Features of the SFH of LGS-3 

Figure E] shows plots of the SFH of LGS-3, iplt, Z), as a function of both t and Z derived from 
the two different photometry sets, ^(t) and ip{Z) are also shown in the ip — t and ijj — Z planes, 
respectively. The projection of i{j{t, Z) on the age-metallicity plane shows the AMR, including the 
metallicity dispersion. 

In Figure El we show ilj{t), the AMR, and the cumulative mass fraction of LGS-3, as a function 
of time with their associated errors. Three vertical dashed lines show the ages of the 10th, 50th, 
and 95th-percentile of the integral of ip{t). Table [3] gives a summary of the integrated and mean 
quantities for the observed area and estimates for the whole galaxy. 

The main features of the SFH of LGS-3 can be observed in Figure [6l The bulk of the star 
formation of LGS-3 occurred at old epochs, with a main peak located ~ 11.5 Gyr ago. The initial 
episode of SF lasted until roughly 9 Gyr ago. Approximately 85% of the stars were formed in this 
initial episode. Since then, the star formation rate has been very low and slowly declining until the 
present. The SFR at the measured main peak is a factor of 15 larger than the mean SFR between 
and 9 Gyr ago (because of limited time resolution, this is a lower limit, see §4.2p . It is worth 
noting that the age at which the rate of the AMR starts to increase more steeply matches up with 
the age at which the SFR decreases to a low level. This may indicate a scenario in which the stars 
were formed at a higher rate at old epochs from fresh, low metallicity gas. When the initial episode 
of star formation was over, the gas reservoir would have been nearly exhausted. As a consequence, 
subsequent star formation would occur at a much lower rate and from an ISM including a high 
fraction of processed gas, hence significantly increasing the metallicity enrichment rate. 

It is interesting to test whether the number of SNe produced in LGS-3 during its lifetime are 
sufficient to produce the chemical enrichment we have derived. Focusing on the main episode, 
1.7 X 10^ M0 of gas were converted into stars from 13.5 to 8.5 Gyr ago. Assuming a minimum 
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progenitor mass for core collapse SNe of 6.5 (jSalaris Cassisjbood i we obtain a total of 2.9 x 10 
SNe produced in the old episode. Accounting for stellar mass loss, and assuming a minimum initial 
mass of 10 Mq for SNe progenitors, then the total number of SNe decreases to 1.6 x lO'^. This 
can be compared to an estimate of the total number of SNe necessary to produce the observe d 
chemical enrichment. Assuming the solar abundance distribution from iGreyesse &: Noeld (119931) . 
that only SNe II are inv olved, and using the iron mass production rate from lNomoto et al.l ( jl99^ ) 
(see lMonelli et al.ll2010d . for more details), we estimate that ~ 80 SNe II are enough to produce an 
increment of metallicity from to 0.0007, which is the upper value obtained for LGS-3 8.5 Gyr ago. 
This estimate was derived assuming the iron mass production of a 13 Mq star (which may vary by 
a factor of 2-3 with different assumptions). Although relatively uncertain, our calculation shows 
that the galaxy needs to retain only roughly 1% of the material released by SNe in order to produce 
the observed metallicity increase. The same calculation can be made for the chemical evolution 
from 8.5 Gyr ago to the present. In this case, the mass of gas converted into stars is 2.4 x 10^ Mq 
and the corresponding numbers of SNe are 4.1 x 10^ and 2.4 x 10^. In this time, the metallicity 
increases from Z = 0.0007 to a current upper value oi Z = 0.003, requiring 290 SNe. This is less 
than 15% of the total, a much larger fraction than calculated for the initial star formation episode. 
This implies that a significantly larger fraction of the newly synthesized material was captured at 
later times at lower SFRs. Note that only SNe II have been considered for this estimates; if SNe la 
were also included, even lower fractions of captured newly synthesized material would result. Note 
also, that much of the latter enrichment could be achieved by SNe la from the initial star formation 
episode. 

The SFH shown in Figure [6] corresponds to the observed area, which covers the center and a 
major fraction of the galaxy. The total mass of stars created in this area over the lifetime of the 
galaxy is 2.0 x 10^ Mq. Assuming an exponential density distribution of scale-length 112 pc (S. 
Hidalgo et al., in prep.), we estimate that the total mass of stars ever formed in LGS-3 is 1.1 x 10^ 



Mq. Adding in the gas mass (| Young &: Loll 199 71 . see 31]) ! then the total baryonic mass of LGS-3 is 
1.5 X 10^ Mq. These values correspond to the best solution with the adopted IMF (x = 1.3 for the 
interval 0.1 <m< 0.5 Mq; x = 2.3 for the interval 0.5 < m < 100 Mq; see Table [2]). If a Salpeter 
IMF with an exponent of 2.35 for the range 0.1 to 100 Mq is used, then the result would be larger 
by a factor of ~ 1.5. 

Figure \7\ shows the observed CMD, the best-fit CMD, and the corresponding residuals. The 
calculated CMD has been built using lAC-star with the solution SFH of LGS-3 as input. Resid- 
uals are given in units of Poisson uncertainties, obtained as {n° — n'^)/y^r}^, where n° and are 
the number of stars in bin i of an uniform grid defined on the observed and calculated CMDs, 
respectively. 
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4.2. Estimating the effects of time resolution on the SFH 

Observational uncertainties as well as the uncertainties inherent to the SFH computation 
procedure result in a smoothing of the derived SFH and a decrease of the age resolution (see 
appendix IA.3P . Having a accurate estimate of this age resolution is fundamental in order to obtain 
reliable information of the actual properties of the galaxy. We have calculated the age resolution 
as a function of time by recovering the SFH of a set of mock stellar populations corresponding to 
very narrow bursts occurring at different times. Five mock populations have been computed, each 
one with a single burst at ages of 1, 4, 8, 11.5, and 13 Gyr. We used lAC-star to create the CMDs 
associated with each mock stellar population. For the mock populations at ages of 11.5 and 13 Gyr, 
five simulations have been done with different random number seeds in order to better characterize 
them. The assumed metallicities and SFRs correspond to those obtained for LGS-3. Observational 
uncertainties were simulated in the CMDs by using the results of the DOLPHOT completeness 
tests. The SFHs of the mock stellar populations were then obtained using the identical procedure 
as for the real data. In all cases, the recovered SFH is well fitted by a Gaussian profile with a 
depending on the burst age. Figure [8] shows the input bursts and the corresponding Gaussian 
profile fits to the solutions, including their peak age and a. From these results we estimate an 
intrinsic age resolution of ^ 1.1 Gyr at the oldest ages which decreases to ~ 0.5 Gyr at an age of 
1 Gyr. Also, a slight shift is observed between the input and solution ages. For the oldest trials, 
the shift is toward younger ages (up to 0.4 Gyr for the oldest trial) and for the youngest trials it 
is toward older ages (about 0.2 Gyr). 

It is particularly important to determine the oldest age at which star formation occurred in 
LGS-3 and the real duration of the first, main star formation episode. Our SFH solution (Figure 
El is well fitted by a Gaussian profile with a peak, iJ.obs = 11.5 Gyr and aobs = 1-2 Gyr in the age 
range 9.5-13.5 Gyr (using the DOLPHOT photometry). The shifts shown in Figure [8] indicate that 
the actual star formation peak would have occurred at about 11.7 Gyr ago. A comparison of the 
solution with our estimated resolution indicates that the initial episode of star formation is just 
barely resolved in our observations. A quadratic subtraction of the cr ~ 1.1 that we have obtained 
for a sharp burst at this age from the derived aobs = 1-2 provides a first estimate of the actual, 
intrinsic dispersion of (Jreai ~ 0.5. 

Given the importance of this measurement, we have pursued a more detailed analysis. We 
have extended the former tests to a set of mock stellar populations with SFHs defined as Gaussian 
profile shaped bursts of different a and centered at 11.5 Gyr. Observational uncertainties were 
simulated in the associated CMDs and the SFHs were recovered following the same process as for 
the galaxy. In all cases, the recovered SFHs are well fitted by Gaussian profiles in the age range 
9.5—13.5 Gyr. Figure [9] shows a quadratic fit to the standard deviation of the solutions, 
a function of the input ones, cr^°'^^. If the actual LGS-3 SFH were well described as a Gaussian 
profile, its a could be inferred by interpolating in this fit, which gives afj^^ = 0.6 Gyr. We will 
adopt this value for the forthcoming discussion. 
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4.3. Comparison with other SFHs derivation methods. 



As for the rest of the LCID galaxies , we h ave obtained the SFH of LGS-3 using the Match 
(jPolphinl I2OO2I ) a nd Cole (ISkillnian et al.ll2003l ) methods in addition to the lAC method. These 



two methods use Girardi et al 



pOOg ) for the stellar evolution library and we have adopted the 



DOLPHOT photometry to be used with them for this comparison (for a descriptio n of the main 



features of these methods and the particulars of their application to LCID galaxies, see lMonelli et al 



201 Obi ). For a proper comparison, we have also obtained the SFH with lAC-star/MinnlAC/IAC- 
pop using the same library and photometry. This also allows us to compare the solutions obtained 
with the I AC method using two different stellar evolution libraries. 

Figure [10] shows the cumulative mass fraction of LCS-3 as obtained with the three tested 
methods. The agreement between the three methods resulting in the same cumulative mass fraction 
(0.8) at 8 Gyr is very good. This strongly supports the conclusion the majority the star formation 
in LCS-3 occurred during this initial episode of star formation. Also, all three methods find star 
formation continuing at much lower rates, until the present. The similar behavior of the three 
solutions allows us to be confident that our scientific conclusions are independent of the SFH 
solution method. In particular, the fact that the galaxy has formed most of the stars after the 
reionization epoch is also confirmed and even reinforced. 

Figure [10] also shows effect of using different stellar models by comparing the lAC-star/MinnlAC/IAC- 
pop results using two different stellar libraries. Here we see that the SFH derived using the 
BaSTI stellar evolu tion library has a slightly faster rise compared to the SFH derived using the 
Girardi et al.l (I2OOOI ) library. However, otherwise, the two are in very good agreement, and the 
main features of the SFH of LGS-3 are found in both derived models at roughly equivalent levels. 

An analysis of the actual resolution provided by the Match and Cole methods, of the kind 
of that done for lAC-star/MinnlAC/IAC-pop in §4.21 would be necessary before going deeper in 
this comparison, which is beyond our objectives here. In the following we will use only the results 
obtained with lAC-star/MinnlAC/IAC-pop and the BaSTI stellar evolution library. 



5. The SFH of LGS-3 compared to other LCID galaxies 



In this section, we will compare th e SFH of LGS-3 wit h that of Phoenix (the other dSph/dIrr 
transition galaxy of the LCID sample. iHidalgo et al.l l2009l ) and with those of Cetus and Tucana 



(the two dSphs of the sample, iMonelli et al 



2010b 



Figure [TT] shows the SFH and the AMR of LGS-3 compared to those of Phoenix (IHidalgo et al 



2OO9I ). The correspondence between the main properties of both galaxies is remarkable. Both 
galaxies show a dominant period of initial star formation at early times, with a result that both 
galaxies had formed 50% of their total stellar mass about 10.5-11.0 Gyr ago, or before a redshift 
of 2. Both galaxies had also formed more than 95% of their total stellar mass by 2 Gyr ago. The 
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AMRs of Phoenix and LGS-3 show similar behaviors from the earhest times, with [M/H] ~ —1.7, 
until reaching 95% of the cumulative mass fraction, with [M/H] ~ —1.2. Note that, in both cases, 
the slope of the AMR increases about ~ 8 — 9 Gyr ago, after the initial episode of star formation. 

The SFHs of both transition type galaxies are similar to those of the two dSphs in the sense 
that an old, main episode exists and that a sharp drop off in the star formation rate occurred ~ 9 
Gyr ago. After this, the transition galaxies. Phoenix and LGS-3, were able to for m stars at low 



steadily decreasing rates while the dSphs completely stopped their star formation. iMonelli et al 



(j2010d ) showed that the first, main episode of star formation occurred earlier in Tucana than in 
Cetus. The main episode of star formation LGS-3 has an age comparable to that of Cetus and is 
also clearly younger than that of Tucana. 

There are significant differences between the AMRs of the transition and the dSph galaxies 
even for t > 9 Gyr. The metal enrichment from the earliest times up to ~ 9 Gyr ago is less than 0.1 
dex in LGS-3 and Ph oenix, while it is more than 0.4 dex in the same period in the dSph galaxies 
(iMonelh et al.1 l2010bl l3). In the case of the transition galaxies, more than 80% of the stars 



were 

formed with roughly the same metallicity. This result is supported by the narrow RGB of LGS-3 
as described in f|2l indicative of a low metallicity spread for the intermediate-to-old age stars. 

The differences between the AMRs of the transition and dSph galaxies is intriguing. Galaxy 
mass is the fundamental parameter which is thought to determine a galaxy's ability to retain gas. 
This might naively imply that the dSphs were initially more massive systems than the transi- 
tion galaxies. However, today, Cetus, Tucana, LGS-3, and Phoenix all have comparable lumi- 
nosities. Tidal strippi ng is thought to be very important in the evolution of dSph galaxies (see 



Peharrubia et al.l l2008l . and references therein). The comparable present day luminosities, com- 
bined with the difference in AMRs, could point to significantly larger halo masses for the isolated 
dSphs which have been lost to tides, or, in other words, a much smaller impact of tides on the 
evolution of the transition galaxies. 



6. Cosmological considerations and the early evolution of LGS-3 

6.1. Background 



Dwarf galaxies are the focus of a major cosmological problem affecting the CDM scenario; the 
number of dark matter sub-halos around Milky Way type galaxies predicted by CDM simulations 
is ari order of magnitude larger than observed ([Kauffmann et al.l Il993l : iKlvpin et al.l Il999l : iMoore 
19991 ). Several explanations have been proposed to solve this problem but most of them use two main 
processes that can dramatically affect the formation and evolution of dwarf-sized halos: heating 
from the ultraviolet (UV) radiation arising from cosmic reionization and internal supernovae (SNe) 
feedback. Both are, in principle, capable of stopping the star formation in a dwarf halo and even to 
fully remove its gas. The detailed information about the SFH at early and intermediate ages, such as 
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we have obtained here, may provide fundamental insight into the problem and help to discriminate 



can be found in Mac Low &: Ferrara (1199 



m02). iKravtsov. Gnedin. fc Klvpin 



(|200J) 



Bullock. Kravtsov. & Weinberg (2001). 


Stoehr et al. 


Ricotti &: Gnedin 


( 


2005 


), 


Strieari et al. 


( 


2008), and 



Sawala et al.l (j2010l ). among others. Here we will give a short summary just to place our forthcoming 
discussion in context. 

Feedback from S Ne can suppress star form ation or even blow out the gas completely from 
the smallest systems. iMac Low &: Ferraral (| 19991 ) presented detailed models showing that the mass 
ejection efficiency is low in galaxies with baryonic mass above ~ 10^ Mq and that only galaxies with 
baryonic mass below ~ 10® Mq could have their gas completely blown away almost independently 
of the SNe mechani c al lu minosity. It is worth mentioning that, according to the assumptions of 
Mac Low Sz Ferraral ( 1999 ). a baryonic mass of ~ 10® Mq would correspond to a total mass of 
6.8 X 10^ Mq, while a baryonic mass of ~ 10^ M0 would correspond to 3.5 x 10^ Mq. 

There is a consensus that UV background heating establishes a characteristic time-dependent 



minimum mass (a filtering mass) for halos that can accrete gas (e.g., lGnedinll200d : IKravtsov et al 
2004 ). The UV background will heat the gas in low mass halos, preventing star formation, and 



eventually photo-evaporate them. As a result, halos not massive enough to have accreted gas before 
the reionization epoch would fail to do so later, unless they accrete mass faster than the rate at 
which the filtering mass increases. The red-shift of the reionization epoch has been estimated from 
polar ization observations of the CMB to be 2 = 10.9 it 1.4 (WMAP 5-year results. iKomatsu et al. 
2OO9I ). However, the presence of the Gunn-Peterson trough in quasars at z ~ 6 shows that the 



universe was not yet fully re- ionized at earlier epochs ( Loeb Sz Barkana 2001 : Becker et ahlboOll ). 



According to models (see e.g.. lMac Low &: Ferraralll999l ). the minimum circular velocity for a 
dwarf halo to accrete and cool gas in order to produce star formation is 30 kms^^, which cor- 
responds to a total mass of ~ 10^ Mq. However, most dwarf galaxies in the Local Group, including 
all the Milky Way satellites and five out of the six of the LCID sample, show circular velocities 
well below that value. Nevertheless, t here a re currently several proposed ways to overcome this 
apparent contradiction. iBullock et al.l (j200ll ) has proposed that dSph and the rest of the smallest 
halos would have been formed during the pre-reionization era. The processes of photoionization 
feedback resulted in the suppr ession of the star f ormat ion below an observable level in 90% of 
them. This was supported bv iRicotti &: GnedinI (|2005l ). who found that if positive feedback is 
considered, a fraction of low mass halos (total mass below 10^ Mq) would have been able to form 
stars before the reionization era. The star formation in these galaxies would have been halted by 
internal mechanisms, like photo-dissociation of H2 or SNe feedback in advance of the reionization 
era. Another possibility is that masses of dark matter halos of the dSphs galaxies would be much 
larger than those measur ed at the optical lim it of the galaxies and larger than the limit required by 
the reionization scenario (IStoehr et al.ll2002l ). In these scenarios, dark matter halos of lower masses 
would have failed to form stars and remain dark. Alternatively, dark matter halos could have been 
much larger in the reionization epoch, with total masses above the mass required for star forma- 



- 15 - 



tion, but they lost a large fraction of their mass afterward due to tidal harassment (jKravtsov et al 
2004 ). It is also possible that a self-shielding mechanism would be present, protecting the gas in 
the inner regions or in r egions denser than some limit from being heated by the UV background 
(jSusa &: Umemural 12004 ) ■ Finally, it has been suggested that inhomogene ous reionization ca n lead 
to large variations in the masses of halos that survive and produce stars (jBusha et al.ll2010l ). 



In summary, it is usefu l to note that the first scenario (IBullock et alj l200ll : iRicotti &: Gnedin 



20051 ) and also the last ones (jSusa &: Umemurall2004l : iBusha et al.ll20ld ) , allow star formation in low 
mass halos (total mass below 10® — 10^ M©), while the second and third are compatible with star 
formation in more massive halos only (total mass above 10^ Mq). Additionally, the first scenario 
would imply that star formation would have been produced before reionization and would have 
been halted afterward in the smallest galaxies. 



Recently, ISawala et al.1 (|20ld ) have presented a set of high-resolution hydrodynamical simula- 
tions of the formation and evolution of isolated dwarf galaxies including the most relevant physical 
effects, namely metal-dependent cooling, star formation, feedback from Type II and la SNe, and 
UV background radiation. Their results are very useful for a direct comparison with observations. 
They study halos with present day total masses between 2.3 x 10*^ Mq and 1.1 x 10^ Mq and reach 
an interesting set of conclusions that include many of the aforementioned ones. First, halos that are 
not massive enough to accrete sufficient gas to form stars before z = 6 lose their gas subsequently 
due to UV background heating. In this sense, reionization sets a lower mass limit of the halo for the 
star formation. Second, in halos that form stars, feedback is the main process driving the evolution 
of the galaxy and regulating its star formation. It alone can blow out all the gas of a dwarf system 
before z = 0, while UV background alone has almost no effect. However, if feedback has previously 
made the gas diffuse and reduced its radiative cooling efficiency, then the UV background has a 
strong effect, producing a sharp cut off in the star formation at the epoch of reionization. Finally, 
self-shielding is effective only in the inner regions of the halos at the more massive end. 

All the aforementioned mechanisms are founded on solid and extensive theoretical modeling. 
However to decide to what extent one or another should be preferred, detailed observational evi- 
dence is necessary. One such observable is whether the star formation in the smallest galaxies was 
halted at reionization. To answer this que stion, a precision of ^ 1 Gyr is required in the determi- 
nation of the SFH at 12.5 — 13 Gyr ago (jRicotti fc GnedinI l2005l ) . which is indeed what we have 
reached in our work. In the following we will discuss the different possibilities in light of our results 
for LGS-3 paying attention first to SNe feedback effects and second to UV background effects. 



6.2. SNe feedback in LGS-3 



The SFH obtained for LGS-3 allows us to make a direct estimate of the importance of feedback 
in LGS-3. To do so, the mechanical luminosity of the SNe produced in the main star formatio n 
episode and the mass of the galaxy can be compared with the results of lMac Low &: Ferraral (|l999l ) . 
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They used models to calculate the conditions for blow-away, partial mass loss, and no mass loss of 
dwarf galaxies as a function of the baryonic mass of the galaxy and the mechanical luminosity, Ly^j , 
of the SNe produced during a central, instantaneous star formation episode. 

The kinetic energy released to the interstellar medium by the SNe produced in the old, main 
star formation episode in LGS-3, can be coi nputed as follows. For a continuous star formation rate 
of 1 Mq yr~^, the STARBURST99 models (jLeitherer et al.lll999l ) predict a mechanical luminosity 
of logL^ = 41.9. This is the equilibrium value reached after a few tens of million years from the 
first star formation. The star formation rate is for a mass range 1 to 100 Mq and is computed 
using a classical IMF with exponent —2.35. For direct comparison, we convert this rate to the same 
stellar mass interval and IMF that we have used and obtain 1.7 Mq yr~^. Scaling these values to 
the SFR obtained for LGS-3, the total mass converted into stars in the initial episode is 1.7 x 10^ 
Mq. Assuming a Gaussian profile of a = 0.6 Gyr, or a duration (FWHM) of 1.4 Gyr, results in an 
average SFR of 1.3 x lO'^ Mq yr^^. Scaling STARBURST99 models to this average SFR results 
in a mechanical luminosity of = 5.3 x 10^^ erg s^^. 

Alternatively, we can compute the mechanical luminosity from core collapse SNe only (note that 
STARBURST99 includes all sources of mechanical luminosity). From ^ the total number of SNe 
formed during the old, main episode was 2.9x10^ or 1.6x10"^ for minimum core collapse progenitor 
masses of 6.5 Mfr^ or 10 M^r^, respecti vely (jSalaris fc Cassisill2005l ). Assuming an energy release per 
SN of 10^^ erg (ILeitherer et al.lll999l l. we obtain a total mechanical luminosity released during the 
old, main episode of = 6.8 x 10^^ erg s~^ or = 3.8 x 10^^ erg s~^, respectively. Since these 
values are compatible with the STARBURST99 estimate, we will adopt 



5.3 X lO^'^ erg s~^ 



Comb ining this with a baryonic mass of 1.5 x 10^ Mq (S ]4.ip . the model results of lMac Low fc Ferrara 
(j 19991 ) shown in their Figure 1, place LGS-3 in the regime of mass loss, but fairly close to the blow- 
away regime. 



However, equation 1 from iMac Low &: Ferraral (| 19991 ) gives a relation between the baryonic 
and dark i natter mass of dwarfs which is an extrapolation to low mass galaxies of the relation 
derived by Persic. Salucci. fc Stel (jl996l ) for more massive galaxies (total masses _above IO^Mq). 
From this relation, the total mass for LGS-3 is 9 x lO"^ Mq. Alternately, ISawala et al.l ( 20ld ) give 
the final {z = 0) stellar, gaseous, and total masses for their model runs. The relations between 
these parameters depend on the model ingredients, but for a baryonic mass like that of LGS-3, 
the total mass is i n the range from 2 x 10^ Mq to 4 x 10^ Mq. Using this estimate instead of the 
Persic et al.l (jl996l ) extrapolation places LGS-3 well within the regime of mass-loss due to feedback 
and far from the blow-away regime. 



6.3. The effects of reionization 



A flat Einstein-de Si tter universe with Hn = 70.5 km s Mpc , ftm = 0.273 (as derived from 
the WMAP 5-year data; komatsu et ahlbood ). has an age ~ 13.7 Gyr. The WMAP 5-year data 
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indicate the epoch of reionization occurred at z ~ 10.9 corresponding to a look-back time of ~ 13.3 
Gyr. Evidence from quasar spectr a indicate the epoch of reionization is over at a redshift of z ~ 6 



or a look-back time of ~ 12.7 Gyr ([Becker et al.ll200ll ). We will focus on the latter since it is more 



conservative from the point of view of the forthcoming discussion. 

Did LGS-3 form most of its stars before or after the epoch of reionization? Using the cumulative 
mass fraction from Figure [6l we find that more than 80% of the stellar mass of LGS-3 has been 
formed later than ~ 12.5 Gyr ago (z ~ 5). Thus, the observations indicate that the majority of the 
stars in LGS-3 were formed after the epoch of reionization. 



We can use synthetic data to better assess this constraint. In Figure \12\ we show the results 
of two synthetic SFHs, one in which all stars are formed between 13.4 and 12.4 Gyr ago (before the 
epoch of reionization) and another where all of the stars are formed between 11.9 and 10.9 Gyr ago. 
Figure [T2] shows that for the model where all of the stars formed before the epoch of reionization is 
clearly inconsistent with the SFH of LGS-3. On the other hand, the model with all stars forming 
after the epoch of reionization is a good match to our SFH for LGS-3, indicating that it is highly 
unlikely that LGS-3 formed a large fraction of its stars before reionization had concluded. 

However, the absolute ages in our comparisons are intrinsically dependent on the stel lar evo- 
lution models which have intrinsic uncertainties (e.g.. lChaboyeiill995l : IChaboyer et al.lll998l ). Mod- 
ern stellar evolution models ha ve addressed many of the shortcomings in earlier generations of 
models (e.g., see discussion in iDotter et al.l 120071 ). and the results of the relative ages derived 
from main sequence stars are in excellent agreement (better than 0.5 Gyr) as demonstrated by 
Marm-Franch et al.l (|2009l ). The accuracy of the absolute ages is also increasi ng. For example, the 
average age of the oldest globular clusters obtained from the age analysis by iMarm-Franch et al 



(|2009l ). using the BaSTI library, is 12.3 Gyr with an age dispersion of ±0.4 Gyr (Marm-Franch, 
private communication). In this regard, the complete inconsistency of our SFH for LGS-3 with the 
model in which all stars are formed before reionization is persuasive. 



Nonetheless, as argued in lMonelli et al.l (|2010d ). the strongest constraints come from relative 
ages. Because all of the LCID galaxies were observed and analyzed in an identical manner we can 
make the following argument. Since we find an age difference of > 1.5 Gyr between the peak of 
star formation in LGS-3 and that of Tucana, and since that difference is larger than the difference 
between the age of the universe and the epoch of reionization, it is highly unlikely that a significant 
fraction of the stars in LGS-3 have been formed before the ep och of reionization. An identical 
conclusion was found for the Cetus dSph (jMonelli et al.l l2010bl ) . In summary. Figure [5] provides 
strong evidence that most, if not all, of the star formation in LGS-3 occurred after the end of 
reionization. 

Although LGS-3 is a low mass galaxy, do theoretical models expect it to form its star before 
reionization? Before putting the former result in the context of current theoretical model pre- 
dictions, some estimates of integrated physical quantities of LG S-3 are necessary. The velocity 
dispersion of the gas in LGS-3 is 8 kms~^ at a radius of 470 pc ( Young Sz Lol 119971 ) which corre- 
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sponds to a total mass within that radius of 2.1 x 10^ Mq. According to IStoehr et al.l ( 20021 ) this 
value should not be considered indicative of the total mass of the ga laxy, but only of the virial mass 
at the measured radius. Following the results of lSawala et al.l (|20ld ). the total mass of LGS-3 would 
be expected to be in the interval ~ 2 X 10^ Mq to ~ 4 X 10^ Mq. Equat ion 16 of lMac Low Ferraral 
gives a relation between circular velocity and total mass in a density distribution given by 

1 /3 

an modified isothermal sphere. Using h = 0.7, it can be written as Vc{r) = 10.5 , where M^j 
is the total mass in units of 10'' Mq and Vc is given in km s^^. Thus, for a total mass in the range 
of 2 X 10® Mq to 4 X 10® Mq, the resulting circular velocity for LGS-3 is between 28 km s ^ and 
36 km s~^. These values are close to or somewhat above the limit of 30 km s~^ below which the 
UV background is expected to prevent any subsequent star formation. 

In summary, it is possible that LGS-3 was massive enough that its star formation was not 
halted by UV background heating and i t was also able to p reserve at least part of its gas against 
feedback from SNe. However, models bv lSawala et al.l (|20ld ) indicate that, for its estimated mass, 
the combination of the the UV background together with SN feedback should have halted the star 
formation at z ~ 6 (or ~ 12.7 Gyr ago). Since this is evidently not the case, a mechanism such 
as self-shielding seems necessary to allow t he galaxy continue forming stars, at least in its inner 
regions, as suggested bv lSawala et al.l (j2010l ). If the total mass of LGS-3 would be of the order of 10® 
Mq or lower (as obtained from the relation used bv lMac Low &: Ferraralll999l ). the circular velocity 
would drop to values of the order of 20 km s~^. For this total mass estimate, our observation that 
most of the star formation in LGS-3 has taken place after z ~ 6 would be difficult to reconcile with 
models. Ho wever, one possibility is that a mechanism like mass loss due to tidal harassment (as 
propo sed by lKravtsov et al.ll2004l ) would be at play. On the other hand, the iMac Low &: Ferrara 
(|l999l ) r elation requir e s a la rge extrapolation for the mass of LGS-3, and this may simply imply 
that the ISawala et al.l (120101 ) mass estimates are to be preferred. 



7. Summary and Conclusions 

We have presented an analysis of the SFH of the transition-type dwarf galaxy LGS-3, based 
on deep HST p hotometry obtained with the ACS. The assumed distance of the galaxy is 652 kpc 
(|Bernardll2009l l. A summary of the results follows. 



• A deep CMD reaching the oldest main-sequence turn-offs with completeness of ~ 0.90 has been 
obtained. Exhaustive crowding tests have been performed in order to properly characterize 
the observational effects. Careful simulations of the observational effects have been conducted, 
taking into account the variable amount of crowding within the galaxy. 

• The SFH for the entire lifetime of the galaxy has been obtained. Tests of self-consistency, 
uniqueness, and stability of the solution have been performed together with tests exploring 
the dependency of the solution on the photometric reduction package, stellar evolution library, 
and the SFH derivation code. 
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• The solution shows that the SFH of LGS-3 is dominated by an old, main episode with maxi- 
mum occurring ~ 11.7 Gyr ago and a duration estimated in 1.4 Gyr (FWHM). Subsequently, 
LGS-3 has continued forming stars until the present at a much lower rate. 

• The total mass of stars produced is 2.0 x 10^ Mq. Of that, 1.7 x 10^ Mq corresponds to the 
old, main star formation episode. The cur rent mass in stars and stellar remnants is 1.1 x 10^ 
Mq. Using a gas mass of 3.8 X 10^ M0 (jYoung fc Lolll997l l the total baryonic mass of the 
galaxy is 1.5 x 10*5 Mq. 

• There is little evidence of chemical enrichment during the first main episode of star formation 
lasting until ~ 9 Gyr ago. After that, the metallicity increased more steeply to a present day 
value of Z ~ 0.0025. The epoch in which the enrichment rate increased is coincident with the 
epoch at which the star formation drops to a low value. This suggests a scenario in which 
LGS-3 formed stars mainly from fresh gas during the first part of its life. After 9 Gyr ago 
much of the gas supply was exhausted and stars formed from a larger fraction of recycled gas. 
This resulted in both a lower star formation rate and an increase of the chemical enrichment 
rate. 

• The difference between the lack of early chemical enrichment in the isolated transition galax- 
ies compared to the significant early chemical enrichment in the isolated dSph galaxies of 
comparable luminosities may be indicative that the dSph galaxies experienced significantly 
more tidal stripping. 

• Most or all the star formation was produced in LGS-3 after the reionization epoch, assumed 
to occur ~ 12.7 Gyr ago. 

• The mass of LGS-3 and the mechanical luminosity from SNe associated with the old, main 
episode (L^ = 5.3 x 10^^ erg s~^) indicate that mass-loss by galactic winds is important, 
but that complete b low-a way was not likely according to the models of SN feedback by 



Mac Low Sz Ferraral (|l999l ). In other words, LGS-3 seems to be massive enough to conserve 



at least a fraction of its gas against SNe feedback. 



• According to models by ISawala et al.l (l20ld ). the total mass of LGS-3 is about 2 x 10^ 
Mf7) to 4 X 10^ MfT), whic h corresponds to a circular velocity of 28 km s~^ to 36 km s^^ 
(jMac Low &: Ferraral Il999l ). These values are close to or somewhat above the limit of 30 km 
below which the UV background radiation would prevent any subsequent star formation. 



According to ISawala et al.l (j20ld ). the combined effects of feedback and UV background 
radiation could have stopped the star formation in a galaxy like LGS-3 after reionization, 
unless some other mechanism l ike self-shielding is present. Indeed, comparison of the SFH 
of LGS-3 with model results by ISawala et al.l (|2010l ) indicate that self-shielding should be at 
play to prevent complete cessation of star formation in LGS-3 at the epoch of reionization. 
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The results presented in this paper allow us to sketch the following scenario for the early 
evolution of LGS-3. The peak of star formation occurred ~ 11.7 Gyr ago, and the star formation 
rate dropped dramatically ~ 9 Gyr ago or earlier although a low star formation activity is observed 
until present. Most or all the star formation occurred in the galaxy after the end of reionization 
(although some star formation before that time cannot be ruled out). Feedback from SNe would 
have created outflows, but would have been insufficient to completely blow-away the gas. The drop 
in star formation at ~ 9 Gyr ago may have been due to the combined results of feedback, UV 
background heating , and gas exhaustion due to star formation. According to models presented by 



Sawala et al.l ((20101), at the epoch of reionization, (in the absence of self-shielding), the combination 



of heating by cosmic UV-background and SNe feedback should have completely stopped the star 
formation in the galaxy. The presence of the star formation activity extending to the present leads 
us to infer that a self-shielding mechanism has been important in the early evolution of LGS-3. 
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A. Confidence of the solution 

As we have mentioned in ^ solutions which are unique, stable, self-consistent, and as inde- 
pendent as possible of the choice of parameters and model functions are desirable. We have done 
many tests to evaluate these properties in the SFHs derived for LCID galaxies. Here we briefly 
discuss these aspects of the calculated SFHs. 



A.l. Uniqueness 



This condition is fulfilled if combinations of simple populations which are within the error bars 
of the adopted solution produce CMDs which are indistinguishable from the best fit CMD. It is also 
required that combinations of simple populations significantly different from the adopted solution 
produce calculated CMDs significantly different fror n the solution one. Simulations done during the 
tests for the lAC-star/MinnlAC/IAC-pop method (|Aparicio fc Hidakol [2OO9I : Iffidalgo et aPbood . 
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and this work), consisting in obtaining the SFH of mock stellar population of known SFH, show 
that the aforementioned conditions are fulfilled in general. 



A. 2. Self-consistency 

By self-consistency we refer to the capability of the method to reproduce synthetic input. When 
a mock population is computed according to a given SFH and analyzed in an identical manner as 
real data, the resulting SFH shou ld be equal to the input on e, within errors and resolution limits. 



Such tests have been executed in lAparicio Hidalgo! (j2009l ) and repeated here using the best fit 



SFH given in Fig. [5] to compute the mock population to be analyzed. For simplicity, we will 
show here only the tests for DOLPHOT photometry but the result is equivalent for DAOPHOT 
photometry. See also tests presented in ^4.2[ 



To minimize the dependency of the test results on the associated statistical fluctuations in- 
herent to the method (present in the observational uncertainties simulation and in the synthetic 
CMD modeling), we have obtained the SFHs for 10 CMDs calculated from the SFH of Figure [5] 
changing only the random number generation seed. An example is shown in Figure [71 Figure [T3l 
shows the average of all of the new SFHs and the AMRs obtained. The input values, i.e., the SFH 
and AMR of LGS-3, are also shown for comparison. The main conclusion of the self-consistency 
test is that the SFH recovering process introduces no significant bias in the location of the main 
peak nor in the AMR, and that the SFH features are properly recovered, although it introduces 
some degree of degradation in time resolution. This has to be taken into account and corrected for 
when quantities like the maximum intensity of a burst or its duration are relevant (see §4.2p . 



A. 3. Stability 

Stability is related to uniqueness, but the question we want to address here is what is the 
maximum excursion of the model functions and external parameters producing a SFH within the 
error bars of the solution. We have tested the effects of changing the fraction of binary stars, the 
IMF, and the photometric zero points. The latter are discussed in ^ The depe ndency of the 



solution on the fraction of binary stars and the IMF has been analyzed in detail in iMonelli et al 



(|2010bl ) and E. Skillman et al. (2011, in prep.) respectively. In short, any fraction of binary stars 
larger than 0.4 gives a SFH within the error bars. In the case of the IMF, for the interval of 
stellar masses 0.5 < m < 100 Mq, the range of values producing solutions within the errors is 
2.2 < X < 2.4, where x is the mass exponent index in ^(m) = Am~^. The solution is independent 
of X for the interval 0.1 < m < 0.5 Mq, although the total stellar mass is not. As an indication, 
the total mass converted into stars in LGS-3 would be increased by a factor ~ 1.5 if the IMF 
exponent changes from x = —1.3 to x = —2.3 for the interval 0.1 <m< 0.5 Mq. The same factor 
of ^ 1.5 exists between the corresponding current masses in existing stars and stellar remnants. 
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As a summary, Table H] shows the range of variation for the binary fraction, IMF, and photometric 
offsets which give solutions within the error bars for the adopted SFH of LGS-3. 

A. 4. Statistical significance 

Since the number of real stars is limited, the observational sampling could be insufficient and 
could introduce undesirable biases in the solution. To check this point we have obtained the SFH 
of five different random subsets of the observed stars. Each one contains 50% of the total, so that 
each star may have been included in several subsamples. Figure [T^ shows the SFH of LGS-3 and 
the average of the SFHs of the five subsamples. For the comparison to be meaningful, the test 
results have been normalized by a factor. Error bars show the dispersion of the five subsample 
SFHs. The agreement is very good and no bias is apparent, which indicates that our solution does 
not suffer from under-sampling effects. 
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Fig. 1.— 



The LGS-3 observed field. Orientation and scale are marked. 
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Fig. 2. — Color-magnitude diagram of LGS-3. In left panel, individual stars are plotted while 
right panel shows star density levels. Star densities bordering the plotted levels run from 8 to 512 
stars dmag^^, evenly spaced by factors of 2. Left axis shows magnitudes corrected from extinction. 
Right axis shows absolute magnitudes. A distance modulus of (m — M)q = 24.07 and an extinction 
Af475w = 0.156, Apsiiw = 0.079 have been used. Bottom lines in the right panel show the 0.25, 
0.50, 0.75 and 0.90 completeness levels. Three isochrones from the BaSTI stellar evolution library 
have been over plotted. 
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2. Error simulation 
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Fig. 3. — Data flow diagram followed in the LCID project to obtain the SFHs of the galaxies. See 
text for details. 
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Fig. 4. — Regions of the CMDs (bundles) used for gridding. Bundles are shown overplotted on the 
observed CMD (left panel) and on one of the synthetic CMD used as an input model for the SFH 
derivation (right panel). Quantities in right panel give the number of bins used to sample each 
bundle. 
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Fig. 5. — Star formation history, 'tp{Z,t), of LGS-3 obtained from DOLPHOT (upper panel) and 
DAOPHOT (lower panel) photometry sets. The SFH as a function of age 1^(1) and metallicity 
il^{Z) are shown projected on the ip — time and ip — metallicity planes, respectively. Dashed lines 
give the error intervals. The age-metallicity relationship is the projection onto the look-back time- 
metallicity plane. 
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Look-back time (Gyr) 



Fig. 6. — (upper panel), metallicity (middle panel), and cumulative mass fraction (lower 
panel) of LGS-3 obtained from DOLPHOT (thick, solid lines) and DAOPHOT (thick, dashed 
lines) photometry sets. Thin lines give the uncertainties. Cyan lines show the final adopted 
solution, which is the mean of DOLPHOT and DAOPHOT ones. Vertical dotted-dashed lines 
indicate the times corresponding to the lO*'^, 50**^ and 95*^^ percentiles of ip{t), i.e. the times for 
which the cumulative fraction of mass converted into stars was 0.1, 0.5 and 0.95 of the current 
one. A redshift scale is given in the upper axis, computed assuming Hq = 70.5 km s~^ Mpc~^, 
$7^, = 0.273, and a flat Universe with Q,a = 1 — ^m- 
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1 2 1 2 1 2 

(F475W- F814W)q 

Fig. 7. — Observed (left panel), calculated (central panel), and residuals (right panel) CMDs. The 
calculated CMD has been built using lAC-star with the solution SFH of LGS-3 as input. The grey 
scale and dot criteria for the first two is the same as for Figure [21 The residuals are in units of 
Poisson uncertainties. 
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[l=^.2 — 

a=0.5^ 




Look-back time (Gyr) 

Fig. 8. — Results of the tests of time resolution. Arrows indicate the burst age of the simulated 
populations. Bursts have a duration of 10^ Mq and their intensities are similar to the SFR of 
LGS-3 at the corresponding age. Arrows heights are not indicative of burst intensity. Curves show 
the Gaussian profile fits to the solutions. The Gaussian profile areas are normalized to the same 
value for easier comparison. The solution peaks (//) and standard deviations (a) are given. 
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Fig. 9. — The standard deviation, cr™"'^ , of Gaussian profile fit solutions obtained for a set of 
Gaussian profile mock input SFHs with standard deviation cr'^'^^, centered at about the age of the 

peak of the main episode of LGS-3 (11.5 Gyr). Observational uncertainties were simulated in the 
associated mock CMDs and the SFHs were recovered following the same procedure as for LGS-3. 
A quadratic fit to the data is shown, together with the a value measured in the real data solution 
("^rec ) the corresponding real value inferred using the fit (o'?^'^). 
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Look-back time 



Fig. 10. — Cumulative mass fraction as a function of look-back time of the solutions obtained 
for LGS-3 with different methods, indicated in the label. All them use the DOLPHOT pho- 
tometry. lAC-B corresponds to the solution adopted and discussed in this paper, obtained with 
lAC-star/MinnlAC/IAC-pop using the BaSTI s tellar evolution libr ary. lAC-G has been obtained 
with lAC-star/MinnlAC/IAC-pop b ut using thelG irardi et al.l (120001 ) stellar evolution library. Cole 
and Match methods make use of the lGirardi et al.. (j200fl) stellar evolution library. 
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Fig. 11. — Compari son between the SFHs and the AMRs of LGS-3 (this paper) and Phoenix 
(jHidalgo et al.l |2009| ). A redshift scale is given on the top axis, computed assuming Hq = 
70.5 km s^^ Mpc~^, = 0.273, and a flat Universe with = 1 — Q^n- Vertical dotted-dashed 
lines indicate the times corresponding to the 10*^, 50*^^ and 95**^ percentiles of ip{t) for each galaxy, 
i.e., the times for which the cumulative fraction of mass converted into stars was 0.1, 0.5 and 0.95 
of the current one. 
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Look-back time (Gyr) 

Fig. 12. — The SFH obtained for two mock galaxies, compared with the actual SFH of LGS-3. One 
of the mock galaxies simulates a system in which most or all the star formation occurred before 
reionization (star formation between 12.4 and 13.4 Gyr ago; magenta, dotted line). The other 

one simulates a system in which all the star formation occurred after reionization (star formation 
between 10.9 and 11.9 Gyr ago; black, dotted line). Magenta and black solid lines show the 
solutions obtained for them while the green solid line correspond to LGS-3. This figure provides 
strong evidence that most, if not all, the star formation in LGS-3 occurred after reionization. 
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Fig. 13. — Consistency test of the SFH of LGS-3. Thick solid hne gives the solution SFH and 
metallicity of the galaxy using DOLPHOT photometry. Thick dashed line shows the recovered 
SFH when the solution SFH is used as mock SFH input. Thin lines show error intervals in top 
panel and real dispersion in bottom panel. 
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Look-back time (Gyr) 

Fig. 14.— Statistical significance of the SFH of LGS-3. Solid lines shows the SFH of LGS-3 and 
its error interval. Thick dashed line shows the average of the SFHs of five subsets of stars of LGS-3 
each one containing a 50% of stars randomly selected. The results of the tests were normalized by 
scaling by a factor of 2. Thin dashed lines show the rms deviations of the five solutions. 
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Table 1. Summary of observations 



Date / UT Exp. Time (s) Filter 



2005- 


■09-12 


16; 


;18; 


;15 


2512 


F475W 


2005- 


■09-12 


16: 


;42; 


;00 


2304 


F814W 


2005- 


■09-14 


16; 


;16; 


;02 


2512 


F475W 


2005- 


■09-14 


16; 


;39; 


;47 


2304 


F814W 


2005- 


■09-15 


14; 


;38; 


;40 


2512 


F475W 


2005- 


■09-15 


15; 


;02; 


;25 


2304 


F814W 


2005- 


■09-16 


13; 


;01; 


;46 


2512 


F475W 


2005- 


■09-16 


13; 


;25; 


;31 


2304 


F814W 


2005- 


■09-17 


11; 


;24; 


;47 


2512 


F475W 


2005- 


■09-17 


11; 


;48; 


;32 


2304 


F814W 


2005- 


■09-17 


14; 


;36; 


;38 


2512 


F475W 


2005- 


■09-17 


15; 


;00; 


;23 


2304 


F814W 
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Table 2. Functions and parameters tested and adopted for the SFH 



Function/Parameter 


Min. value 


Max. value 


Step 


Adopted 


Stellar evolution library 


BaSTI 


Girardi 




BaSTI 


Binary fraction 


0.0 


1.0 


0.2 


0.4 


IMF (0.1 < m < 0.5 Mq)'' 


0.0 


4.0 


0.1 


1.3 


IMF (0.5 < m < 100 Mq) 


0.0 


5.0 


0.1 


2.3 




-0.21 


0.21 


0.105 


0.13/0.08 




-0.15 


0.15 


0.075 


0.11/0.11 


CMD gridding<^'<* 


0.01 X 0.2 


0.01 X 0.2 






Age sampling (< 1 Gyr)'' 


0.2 Gyr 


0.5 Gyr 






Age sampling (> 1 Gyr) 


1.0 Gyr 


1.0 Gyr 






Metallicity sampling (Z < 0.001)'* 


0.0002 


0.0002 






Metallicity sampling (Z > 0.001) 


0.0005 


0.0005 







^Values are the exponent x in the IMF defined as 4)(m) = A m ^, where m is the 
mass and A is a normalization constant. 

Zero point offsets are relative to ISirianni et al. I 1I2OO5I') . The adopted values were 
selected by interpolation. Results for DOLPHOT/DAOPHOT are given in the last 
column. 

'^Values are the box size (color X magnitude) for the CMD sampling. 

Several bin sets were tested using the same bin size but applying different offsets to 
the starting point of the bin set. 
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Table 3. Summary of integrated and mean quantities of the LGS-3 SFH 



Area (kpc^) 


Mass formed(10^ M©) 


<[M/H]> 


< V(t) > (10-3 M0 yr-l) 


0.42 


1.26 ±0.04 


-1.7±0.1 


0.09 ± 0.03 


Total Galaxy 


1.95 ±0.04 


-1.7±0.1 


0.14 ±0.03 
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Table 4. Model functions and external parameters best fit range 



Punction/Pajrameter 


Min./Max. value 


Binary fraction 


0.4/0.8 


IMF (0.1 < m < 0.5 M©)''- 


-4.0/0.0 


IMF (0.5 < m < 100 Mq) 


-3.0/-2.0 




-0.04/0.02 




-0.06/0.04 



^Values are the exponent x in the IMF defined 
as 4>{m) = A , where m is the mass and A is 
a normalization constant. 



